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ABSTRACT 

We present a study of the star cluster population in the starburst irregular 
galaxy NGC 4449 based on B, V, I, and Ha images taken with the Advanced 
Camera for Surveys on the Hubble Space Telescope. We derive the cluster prop- 
erties such as size, ellipticity, and total magnitudes. Cluster ages and masses are 
derived fitting the observed spectral energy distributions with different popula- 
tion synthesis models. Our analysis is strongly affected by the age-metallicity 
degeneracy; however, if we assume a metallicity of ~l/4 solar, as derived from 
spectroscopy of HII regions, we find that the clusters have ages distributed quite 
continuously over a Hubble time, and they have masses from ~ 10 3 M o up to 
~2x 1O 6 M , assuming a Salpeters' IMF down to 0.1 M Q . Young clusters are 
preferentially located in regions of young star formation, while old clusters are 
distributed over the whole NGC 4449 field of view, like the old stars (although 
we notice that some old clusters follow linear structures, possibly a reflection of 
past satellite accretion). The high SF activity in NGC 4449 is confirmed by its 
specific frequency of young massive clusters, higher than the average value found 
in nearby spirals and in the LMC (but lower than in other starburst dwarfs such 
as NGC 1705 and NGC 1569), and by the flat slope of the cluster luminosity 
function (dN(Ly) oc L^^dL for clusters younger than 1 Gyr). We use the upper 
envelope of the cluster log(mass) versus log(age) distribution to quantify cluster 
disruption, and do not find evidence for the high (90%) long-term infant mor- 
tality found by some studies. For the red clusters, we find correlations between 
size, ellipticity, luminosity and mass: brighter and more massive clusters tend to 
be more compact, and brighter clusters tend to be also more elliptical. 

Subject headings: galaxies: dwarf — galaxies: individual (NGC 4449) — galaxies: 
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irregular — galaxies: star clusters: general — galaxies: starburst 
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Introduction 



Star clusters are present in all types of galaxies, from the earliest to the latest 
morphological types, from the quietest to the most turbulent ones, such as merging and 
starburst systems. They are suggested to be the birth site of many (possibly most) stars, 
and provide an excellent tool to study the parent galaxy evolution through a long range of 
epochs. For these reasons, the interest in star cluster studies never fades. While Galactic 
star clusters were already studied a century ago, extragalactic ones have required the 
superb spatial resolution of the Hubble Space Telescope (HST) to receive the deserved 
attention. Nowadays people use HST to (at least partially) resolve both old globulars 
(GCs) and young and intermediate-age clusters up to distances of a few Mpc. Thanks 
to HST, special interest has been focussed in the last two decades to the young massive 
cluster s typical of star f ormin g galaxies, inc l uding t he so-called S u per Star Clusters - SSC s 



O'Connell et al. 



Sabbi et al. 



(J2QQ3); 



(1994); 



Bil 



ett et al. 



Annibali et al 



20091 ): 



(2002 



Larson 



(|2Q09j); 



Larsenl (I2004h: 



Whitmore et al. 



Mora et al. 



(120051 ): 



f j2009h ). These 



objects have masses typical of globular clusters and may thus be the young counterpart 
of the classical fossil records of the earliest localized star formation activity. If globular 
clusters could still be formed today in some environments, their study would provide a 
direct insight into the conditions that were present in the early days of galaxy formation, 
when the globular clusters that we see today in the halos were formed. 

In this paper we present new data on the star clusters in the starburst irregular galaxy 
NGC 4449. NGC 4449 is known to h ost a young (~ 6-15 Myr) centr al super star cluste r 



(IBoker et al. 



2001 



Gelatt et al. 



20011 ). numerous star clusters 



60, 



Gelatt et al. 



( 1200 lh ) 



1 Based on observations with the NASA/ESA Hubble Space Telescope obtained at the 
Space Telescope Science Institute which is operated by AURA Inc. for NASA under contract 
NAS5-26555. 
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with ages up to 1 Gyr and possibly old er, and at least 13 embedded massive clusters 



detected in the radio with ages <5 Myr (IReines et al. 



2008). Our analysis is based on our 



HST data acquired with the Wide Field Channel (WFC) of the Advanced Camera for 
Surveys (ACS). 



The Magellanic irregular galaxy NGC 4449 (a 20 oo = 1 2 fe 28 m ll s .9 5 



I = 136.84 and b = 72.4) at a distance of 3.82 ± 0.18 Mpc flAnnibali et al. 



2ono 



+ 44°05 40 



fl2008f ). hereafter 



A08) is one of the best studied and spectacular nearby starbursts. Its size and metallicity 
are similar to those of the Large M agellanic Cloud (LMC), while its integrated magnitude 



M, 



13 



18.2 flHunter et al 



fll999f ). A08) is « 1.4 times brighter than the LMC's. Also 



its SFR, averaged over a Hubble time, is ~ 0. 2 M^ yr" 1 , twice tha t of the LMC, while 



its current SFR is as high as ~ 1.5 M yr 1 (IThronson et al. 



19871 ). It is indeed one of 



th e most luminous and a ctive irregular galaxies, as confirmed also by the recent study 



by 



McQuinn et al. 



(120 lOh on the star formation histories (SFHs) of 18 nearby starburst 



galaxies, of which NGC 4449 turns out to be the strongest one. 

Data acquisition and reduction are summarized in Section 121 while the cluster analysis 
and the artificial cluster experiments are described in Section |3j The results (luminosity 
functions, ages and masses, cluster classification, correlations between cluster properties) 
are presented in Section HJ In Section [5] we investigate cluster disruption. The final results 
and discussion are presented in Section [6j 



2. Observations and data reduction 

Here we briefly recall the adopted observing strategy and main steps involved in the 
data reduction process, which were fully described in A08. 

The observations were performed in November 2005 with the ACS/WFC using 
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the F435W (B), F555W (V), and F814W (I) broad-band filters, and the F658N (Ha) 
narrow-band filter (GO program 10585, PI Aloisi). We had two different pointings in a 
rectangular shape along the major axis of the galaxy. Each pointing was organized with a 
4 - exposure half + integer pixel dither pattern with the following offsets in arcseconds (0, 
0) for exposure 1, (0.12, 0.08) for exposure 2, (0.25, 2.98) for exposure 3, and (0.37, 3.07) 
for exposure 4. This dither pattern is suitable to remove cosmic rays and hot/bad pixels, 
fill the gap between the two CCDs of the ACS/WFC, and improve the PSF sampling. Four 
exposures of ~ 900 s, 600 s, 500 s and 90 s for each of the two fields were acquired for each 
of the B, V, I and Ha filters, respectively. 

For each filter, the eight frames (4 dithered exposures for each of the 2 fields), 
calibrated through the most up-to-date version of the ACS calibration pipeline (CALACS), 
were co-added into a sin gle mosaicked image using the software package MULTIDRIZZLE 



( IKoekemoer et al.ll2002l ). The MULTIDRIZZLE procedure also corrects the ACS images for 
geometric distortion and provides removal of cosmic rays and bad pixels. The total field of 
view of the resampled mosaicked image is ~ (380 x 200) arcsec 2 with a pixel size of 0.035 
"(0.7 times the original ACS/WFC pixel size). 

The total integration times are ~ 3600 s, 2400 s, 2000 s and 360 s for the B, V, I and 
Ha images, respectively. Only in a small region of overlap between the two pointings (~ 
(30 x 200) arcsec 2 ) the integration times are twice as those listed above. 

The source detection was performed with the DAOFIND task within IRAlj 2 ! above 3 
times the background standard deviation, independently in the three ban ds. PSF-fitting 



photometry was performed with the DAOPHOT package ([Stetson 



19871 ) in the IRAF 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are oper- 
ated by AURA, Inc., under cooperative agreement with the National Science Foundation 
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environment. The PSF was computed from the most isolated stars in the frame, and 
was modeled with an analytic Moffat function plus additive corrections derived from the 
residuals of the fit to the PSF stars. Then we ran ALLSTAR to iteratively fit the PSF 
model to groups of stars, providing coordinates and magnitudes of the individual stars, 
along with magnitude errors, y 2 , and sharpness. In particular, the sharpness parameter 
provides a measure of the intrinsic size of the object with respect to the PSF. Sources with 
profile widths similar to the PSF's have sharpness of about zero, cosmic rays have sharpness 
less than zero, and extended sources have sharpness larger than zero. The PSF FWHM is 
< 0.1", corresponding to ~ 2 pc at NGC 4449' s distance. 

The instrumental magnit udes were transformed into the HST VEGAMAG system 



following 



Sirianni et al. 



(120051 ). Aperture corrections were calculated from isolated stars 



selected in our images. Cor rections for imperfect charge transfer efficiency (CTE), applied 



following the formulation of 



Riess fc Mack 



(120041 ). are negligible for the brightest stars, but 



can be as high as ~ 0.1 mag for the faintest stars. 

The B,V and I catalogues were cross-correlated with the requirement of a spatial offset 
smaller than 1 pixel between the positions of the stars in the different frames. This led to 
299,115 objects having a measured magnitude in both B and V, 402,136 objects in V and I, 
and 213,187 objects photometred in all the three bands simultaneously. 
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3. Cluster analysis 
3.1. Selection 

Cluster candidates were selected from the DAOPHOT photometric catalog adopting 
a magnitude cut of m F si4w ^ 22 [] and requiring sharpnessoAo >0.5 in all three F435W, 
F555W, and F814W filters. After this preliminary selection, we visually inspected the 
objects in all the three images looking for centrally concentrated clusters. At a visual 
inspection, some of the objects turned out to be blends of two or more point-like sources, 
while others appeared as bona fide background galaxies, with irregular fuzzy structures, 
spiral arms, or disky features. These objects were removed from the cluster candidate 
lists. We ended up with a final list of 81 candidates, whose coordinates are provided in 
Cols. 2 and 3 of Table 1. Their location in the NGC 4449 field of view is shown in Fig. [TJ 
Background elliptical galaxies may still be present in the candidate list, but the fact that 
the candidates do not appear uniformly distributed over our field of view suggests that the 
majority of them are indeed clusters belonging to NGC 4449: in fact, most of them are 
either concentrated in the galaxy central bar, or close to the S-shaped filaments delineated 
by stars younger that 10 Myr and by the H Q emission (Figs. 19 and 20 in A08). 

The clusters are individually shown in Fig. [2J We notice that a significant number 
of them are partially resolved into stars in our images, which further confirms that the 
contamination from background sources is minimal. We also show in Fig. [3] color-composite 
images of the clusters associated with Ha emission, exhibiting a variety of ionized gas 
morphologies. Our F658N images are quite shallow, thus it is possible that we are missing 

3 At fainter magnitudes, contamination by NGC 4449 stars, background galaxies, and 
blends becomes severe. Also, our completeness for cluster detection drops significantly, as 
quantified in Section I3T41 
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clusters with very faint ionized gas emission. 



Gelatt et al. 



( 1200 ll ) performed a study of the star cluster population in NGC 4449 
using WFPC2 images. Their data cover a field of view of ~ (240 x 200) arcsec 2 , smaller 
than ours ~ (380 x 200) arcsec 2 (see Fig. QJ. Besides the SSC at the center of NGC 4449, 
they identified 60 candidate compact star clusters, 37 of which are in common with us. We 



Gelatt et al 



(2001) 



visually inspected in our images the objects classified as clusters by 
but not by us. Some of them appear extended in F435W or F555W, but are resolved into 
two or more stars in F814W. Others appear more similar to stellar associations (but we 
can not distinguish between "open clusters" or chance superpositions of unaffiliated stars) 
rather than t o com pact clusters. We show in Fig. 0] two examples of objects classified by 



Gelatt et al. 



(120011 ) as clusters, but considered by us to be "stell ar associations". Fo r the 



Gelatt et al 



(120011). Our 



37 clusters in common, we provide in Table 1 the ID assigned by 
ACS images are saturated in the center of the galaxy, preventing a study of the SSC in 
NGC 4449's center. 



Table 1. Cluster properties 



ID a R.A. Dec R e b e c m F435w d m F555W d m F814w d M v c Age f Mass f My (10 Myr)s 

J2000 J2000 "(pc) Myr 1O 5 M 



ssc 


12 


28 


11 


13 


44 


05 


37 


42 










1 


12 


28 


27 


04 


44 


07 


00 


54 





17 


(3 


2) 


2 


12 


28 


19 


31 


44 


07 


29 


72 





21 


(3 


8) 


3 


12 


28 


16 


44 


44 


07 


29 


49 





49 


(9 


0) 


4 


12 


28 


17 


88 


44 


06 


55 


41 





21 


(3 


9) 


5 


12 


28 


16 


62 


44 


07 


05 


62 





12 


(2 


1) 


6 


12 


28 


16 


87 


44 


06 


55 


72 





33 


(6 


1) 


7 


12 


28 


17 


57 


44 


06 


45 


33 





31 


(5 


8) 


8* (59) 


12 


28 


18 


79 


44 


06 


23 


08 





20 


(3 


7) 


9* 


12 


28 


14 


56 


44 


07 


09 


72 





04 


(0 


7) 


10t* (51)* 


12 


28 


15 


89 


44 


06 


43 


54 





07 


(1 


4) 


11* (58) 


12 


28 


17 


54 


44 


06 


10 


45 





32 


(6 


°) 


12t*(52)* 


12 


28 


16 


05 


44 


06 


29 


64 





08 


(1 


5) 


13* (49)* 


12 


28 


14 


90 


44 


06 


41 


39 





09 


(1 


7) 


14 


12 


28 


11 


65 


44 


07 


20 


74 





08 


(1 


5) 


15* (54) 


12 


28 


16 


69 


44 


06 


12 


18 





20 


(3 


7) 


16 


12 


28 


14 


84 


44 


06 


35 


16 





26 


(4 


8) 


17* (57) 


12 


28 


17 


51 


44 


05 


57 


78 





17 


(3 


2) 


18* (47) 


12 


28 


13 


99 


44 


06 


43 


21 





16 


(2 


9) 


19* (61) 


12 


28 


19 


19 


44 


05 


32 


00 





29 


(5 


3) 


20* (60) 


12 


28 


18 


82 


44 


05 


19 


61 





48 


(8 


8) 


21* (53) 


12 


28 


16 


58 


44 


05 


36 


64 





19 


(3 


5) 


22* (42) 


12 


28 


12 


99 


44 


06 


23 


42 





19 


(3 


5) 






07 


21 


24±0 


01 


20 


49±0 


01 


19 


41±0 


01 





08 


22 


24±0 


03 


21 


52±0 


02 


20 


51±0 


03 





04 


22 


42±0 


08 


21 


63±0 


07 


20 


56±0 


09 





05 


20 


90±0 


01 


20 


52±0 


01 


19 


81±0 


02 





11 


21 


26±0 


01 


20 


97±0 


01 


20 


22±0 


02 





03 


18 


87±0 


01 


18 


59±0 


01 


18 


01±0 


02 





09 


20 


95±0 


04 


20 


19±0 


04 


19 


00±0 


03 





11 


20 


22±0 


02 


19 


48±0 


01 


18 


43±0 


02 





15 


19 


02±0 


01 


18 


92±0 


01 


18 


64±0 


01 





15 


21 


54±0 


01 


20 


22±0 


01 


20 


84±0 


02 





08 


20 


41±0 


03 


20 


21±0 


04 


19 


69±0 


07 





13 


20 


22±0 


01 


18 


55±0 


01 


19 


25±0 


01 





36 


20 


51±0 


01 


20 


63±0 


02 


20 


92±0 


07 





05 


20 


59±0 


01 


19 


87±0 


01 


18 


87±0 


01 





04 


20 


25±0 


02 


20 


01±0 


02 


19 


39±0 


04 





09 


21 


15±0 


04 


20 


85±0 


05 


20 


15±0 


09 





10 


19 


89±0 


01 


19 


70±0 


01 


19 


27±0 


03 





27 


18 


15±0 


01 


17 


83±0 


01 


17 


19±0 


01 





05 


21 


58±0 


02 


20 


73±0 


02 


19 


50±0 


02 





03 


20 


53±0 


01 


19 


75±0 


01 


18 


62±0 


01 





02 


20 


33±0 


02 


19 


59±0 


01 


18 


54±0 


02 





12 


20 


32±0 


02 


20 


07±0 


02 


19 


57±0 


04 



-7.56±0 


15 


1441^4 


„ .„+0 01 


,-.,-,+0 15 


-6.53±0 


15 


+ 285 

36011^5 


^ ~„-i-o 


-11.211°;^ 


-6.42±0 


17 


70351^6 


_ __i_n 03 

0-90lo.o3 


-ll.73l^;i( 


-7.50±0 


15 


6O9I45 


„ — +0 02 

0.37t£o2 


-10.661°;^ 


-7.06±0 


15 


444^| 


0.21+^;^ 


-9.98^ 


-9.42±0 


15 


295_ 10 


1 QC+0.01 

L36 -o.oi 


-12.03_ (hl5 


-7.87±0 


15 


1510lg 4 


„ ^4-0.01 

°- 76 -0.09 


-1-1 in + 0.24 

-11.49+^ 1B 


-8.57±0 


15 


518llggo 


r -1 ri4-0 75 


_j_0.19 

— id.oy_ 22 


-9.08±0 


15 


83 ^28 


_ _ „_i_o 07 

- 50 -oio 


- 10 -87la22 


-7.70±0 


15 








-7.80±0 


15 


259^1 


0.30±°;gi 


-10.35±°;H 


-9.37±0 


15 








-7.33±0 


15 




02+ - 01 


— 7 OR" 1 " 1 ' 83 
' - UD -0.18 


-8.18±0 


15 


3240l|^ 




1 't'7+0. 18 
lz - ' ' -0.18 


-8.01±0 


15 


296±\ 


0.38±°Z 


-10.6210;" 


-7.18±0 


16 


389±°° 


20+ - 02 


-10.00i° ;i 9 


-8.30±0 


15 


251±? 6 


45+ - 01 
u - 4o -0.03 


-10.83±g;i« 


-10.19±0.15 


330^23 
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-7.34±0.15 


13197 1524 
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1450lf 
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Table 1 — Continued 



ID a R.A. Dec R e b e c m F4 , S5 w d m F55BW d m F8 uw d M v c Age f Mass f M v (10 Myr)s 

J2000 J2000 "(pc) Myr 10 5 Mq 



23* (43) 


12 


28 


13 


03 


44 


06 


12 


01 





12 


(2 


1) 


2d 


12 


28 


07 


32 


44 


07 


21 


54 





30 


(5 


5) 


25 


12 


28 


14 


61 


44 


05 


42 


47 





31 


(5 


7) 


26 


12 


28 


13 


23 


44 


06 


00 


10 





21 


(3 


9) 


27 


12 


28 


07 


72 


44 


07 


13 


42 





10 


(1 


8) 


28 


12 


28 


12 


85 


44 


06 


04 


01 





27 


(5 


1) 


29 


12 


28 


14 


45 


44 


05 


41 


87 





12 


(2 


2) 


30 


12 


28 


14 


41 


44 


05 


37 


08 





15 


(2 


7) 


31 


12 


28 


14 


85 


44 


05 


27 


51 





07 


(1 


3) 


32 


12 


28 


12 


81 


44 


05 


50 


06 





38 


(7 


1) 


33 


12 


28 


10 


78 


44 


06 


15 


48 





35 


(6 


5) 


34* (26) 


12 


28 


10 


32 


44 


06 


21 


37 





13 


(2 


5) 


35* (40) 


12 


28 


12 


37 


44 


05 


43 


21 





12 


(2 


3) 


36* (36) 


12 


28 


12 


01 


44 


05 


47 


63 





19 


(3 


5) 


37 


12 


28 


12 


65 


44 


05 


35 


14 





30 


(5 


5) 


38* (37) 


12 


28 


12 


11 


44 


05 


39 


25 





19 


(3 


5) 


39* (48) 


12 


28 


14 


60 


44 


05 


00 


43 





09 


(1 


7) 


40* (30) 


12 


28 


11 


12 


44 


05 


45 


56 





09 


(1 


7) 


41 


12 


28 


12 


25 


44 


05 


29 


89 





24 


(4 


5) 


42* 


12 


28 


11 


30 


44 


05 


42 


41 





21 


(3 


8) 


43*(31)* 


12 


28 


11 


22 


44 


05 


38 


06 





11 


(2 


0) 


44 


12 


28 


11 


74 


44 


05 


28 


03 





36 


(6 


7) 


45* (38) 


12 


28 


12 


15 


44 


05 


18 


15 





13 


(2 


4) 






20 


18 


25±0 


02 


18 


21±0 


02 


17 


71±0 


02 





09 


21 


23±0 


01 


20 


42±0 


01 


19 


34±0 


02 





18 


20 


49±0 


03 


20 


22±0 


04 


19 


60±0 


06 





27 


18 


40±0 


03 


18 


32±0 


03 


18 


04±0 


05 





07 


20 


75±0 


01 


19 


93±0 


01 


18 


81±0 


01 





13 


20 


46±0 


03 


20 


16±0 


03 


19 


49±0 


04 
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20 


82±0 


02 


20 


55±0 


02 


19 


71±0 


03 





08 


21 


05±0 


02 


20 


74±0 


02 


20 


13±0 


06 





10 


21 


70±0 


02 


21 


41±0 


02 


20 


76±0 


04 





10 


20 


66±0 


27 


20 


31±0 


20 


19 


68±0 


20 





10 


20 


67±0 


03 


20 


39±0 


05 


19 


86±0 


09 





24 


19 


40±0 


01 


18 


61±0 


01 


17 


51±0 
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20 


19 


14±0 


03 


18 


39±0 


02 


17 


18±0 
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34 


18 


63±0 


03 
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43±0 
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17 


61±0 
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11 


21 
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05 


20 


74±0 


05 


19 
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06 





23 
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17 
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10 


20 


24±0 


01 


19 


43±0 


01 


18 


27±0 


01 





16 


19 


30±0 


03 


18 


67±0 


02 


17 


59±0 


02 





14 


20 


75±0 


03 


20 


57±0 


04 


20 


13±0 


08 





10 


20 


02±0 


10 


19 


77±0 


10 


19 


09±0 


12 





41 


17 


43±0 


01 


17 


25±0 


01 


16 


93±0 


01 





10 


20 


50±0 


26 


20 


16±0 


18 


19 


44±0 


14 





17 


20 
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02 


20 
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02 
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Table 1 — Continued 



ID a R.A. Dec R e b e c m F4 , 35w d m F555 w d rn F8 i4.w d M v c Age f Mass* M v (10 Myr)s 

J2000 J2000 "(pc) Myr W>M Gl 



46 


12 


28 


10 


24 


44 


05 


42 


57 





18 


(3 


4) 


47 


12 


28 


12 


21 


44 


05 


12 


60 





20 


(3 


6) 


48 (27) 


12 


28 


10 


35 


44 


05 


35 


79 





19 


(3 


5) 


49 1 (16) 


12 


28 


07 


07 


44 


06 


17 


82 





12 


(2 


2) 


50 


12 


28 


09 


21 


44 


05 


40 


42 





04 


(0 


7) 


51 T * 


12 


28 


10 


04 


44 


05 


24 


95 





16 


(3 


0) 


52* (14) 


12 


28 


06 


64 


44 


06 


07 


78 





13 


(2 


4) 


53* (34) 


12 


28 


11 


46 


44 


05 


00 


57 





17 


(3 


1) 


54 


12 


28 


08 


82 


44 


05 


33 


34 





12 


(2 


2) 


55* (23)* 


12 


28 


09 


74 


44 


05 


19 


71 





07 


(1 


3) 


56 1 (22)* 


12 


28 


09 


65 


44 


05 


16 


27 





04 


(0 


7) 


57* (24) 


12 


28 


09 


90 


44 


05 


09 


83 





09 


(1 


6) 


58* (13) 


12 


28 


06 


44 


44 


05 


55 


12 





08 


(1 


4) 


59 


12 


28 


07 


95 


44 


05 


23 


48 





14 


(2 


6) 


60 


12 


28 


09 


59 


44 


04 


58 


43 





25 


(4 


6) 


61 


12 


28 


08 


17 


44 


05 


11 


01 





13 


(2 


4) 


62 


12 


28 


07 


64 


44 


05 


12 


58 





10 


(1 


9) 


63* (17) 


12 


28 


07 


40 


44 


05 


15 


07 





19 


(3 


6) 


64 


12 


28 


10 


24 


44 


04 


36 


28 





05 


(0 


9) 


65* (19) 


12 


28 


08 


82 


44 


01 


52 


01 





26 


(4 


9) 


66 


12 


28 


07 


86 


44 


05 


03 


49 





38 


(7 


1) 


67 


12 


28 


09 


34 


44 


04 


38 


99 





42 


(7 


8) 


68* (11) 


12 


28 


06 


24 


44 


05 


15 


98 





16 


(2 


9) 






06 


20 


65±0 


15 


19 


92±0 


10 


18 


79±0 


08 





17 


21 


22±0 


04 


20 


91±0 


04 


20 


30±0 


08 





04 


18 


68±0 


03 


18 


27±0 


03 


17 


43±0 


02 





18 


20 


47±0 


01 


20 


16±0 


01 


20 


48±0 


02 





27 


20 


99±0 


01 


20 


83±0 


01 


20 


45±0 


01 





09 


20 


46±0 


01 


19 


60±0 


01 


19 


72±0 


01 





13 


19 


95±0 


01 


19 


15±0 


01 


18 


01±0 


01 





11 


20 


00±0 


01 


19 


78±0 


01 


19 


33±0 


03 





19 


21 


53±0 


01 


21 


34±0 


02 


20 


83±0 


03 





33 


19 


21±0 


02 


19 


29±0 


02 


19 


42±0 


05 





26 


18 


09±0 


01 


17 


91±0 


01 


18 


02±0 


01 





05 


18 


38±0 


01 


18 


51±0 


02 


18 


68±0 


04 





12 


19 


65±0 


01 


18 


87±0 


01 


17 


76±0 


01 





11 


21 


23±0 


03 


21 


00±0 


03 


20 


46±0 


06 





14 


18 


97±0 


01 


18 


94±0 


01 


18 


67±0 


02 





22 


20 


67±0 


02 


20 


44±0 


02 


19 


98±0 


05 





13 


22 


70±0 


04 


22 


25±0 


09 


21 


20±0 


09 





09 


20 


97±0 


03 


20 


61±0 


03 


19 


97±0 


06 





26 


21 
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01 


21 


62±0 


01 
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01 





16 


19 
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02 


19 
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02 


19 
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03 





15 


19 
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02 


19 
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02 


18 


73±0 


04 





11 


19 
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01 


19 
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02 


18 
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15 


»±i 
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15 


253+i° 


23+ - 01 
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17 


1265+.JJ 7 
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03 '-167 
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15 
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15 
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15 
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Table 1 — Continued 



ID a R.A. Dec 

J2000 J2000 



Re b 
'(PC) 



mp 555W 



M v c 



Age f 
Myr 



Mass f M v (10 Myr) 
10 5 M Q 



AO 


1 


oo 


177 




i i 


U 1 


01 


^■7 


u 


O 1 


1 1 
. i 


A \ 
J 


U 


Uu 


■71 
L 1 


£1*3 4-17 


177 
U t 


•7n 


Qc;_i_n 
yoztu 


UD 


"1 (7 

ly 


(71 -i-Cl 


(77 


7 


1 n+n 
1U±U 


ID 


1 A OQ^^- 


U.ou_Q 


ru 


12 


28 


n -1 


22 


41 


05 


2 / 


95 





07 


(1 


2) 





06 


21 


97i0 


01 


21 


68i0 


01 


21 


12±0 


02 


-6 


33i0 


15 


2941^ 


0.081q 


71 


12 


28 


03 


09 


11 


05 


31 


92 





32 


(6 


0) 





04 


20 


23±0 


02 


20 


01±0 


03 


19 


59±0 


06 


-8 


00±0 


15 


2511L 


0.331q 


72* (8) 


12 


28 


04 


69 


11 


04 


53 


12 





12 


(2 


2) 





01 


21 


17±0 


02 


20 


36±0 


02 


19 


21±0 


02 


-7 


70±0 


15 


10560:^241 


4.02lg 


73* (15) 


12 


28 


06 


87 


11 


04 


20 


51 





07 


(1 


4) 





12 


19 


49±0 


01 


19 


62±0 


01 


19 


80±0 


02 


-8 


34±0 


15 




0.04lg 


74* (3) 


12 


28 


02 


54 


41 


05 


18 


39 





23 


(4 


3) 





13 


21 


57±0 


06 


21 


40±0 


08 


21 


13±0 


22 


-6 


60±0 


17 




0.06lg 


75* (5) 


12 


28 


02 


61 


41 


04 


51 


51 





26 


(4 


7) 





11 


21 


29±0 


01 


20 


62±0 


01 


19 


58±0 


02 


-7 


43±0 


15 


14271} 


0.47lo 


76* (7) 


12 


28 


03 


88 


41 


01 


15 


26 





17 


(3 


2) 





03 


20 


51±0 


01 


19 


70±0 


01 


18 


52±0 


01 


-8 


36±0 


15 


12000_ 2575 


8.19_! 


77 


12 


27 


57 


53 


44 


05 


28 


16 





10 


(1 


9) 





24 


19 


29±0 


01 


18 


51±0 


01 


17 


42±0 


01 


-9 


55±0 


15 


72441™ 


16.601 


78 1 "* 


12 


27 


56 


63 


41 


05 


29 


32 





15 


(2 


8) 





20 


21 


64±0 


01 


20 


14±0 


01 


20 


73±0 


03 


-7 


78±0 


15 






79 


12 


27 


59 


87 


14 


04 


43 


34 





10 


(1 


9) 





27 


18 


99±0 


01 


18 


19±0 


01 


17 


09±0 


01 


-9 


87±0 


15 


8130lgJ 


24.191 


80 


12 


27 


57 


66 


41 


04 


41 


91 





15 


(2 


8) 





01 


20 


77±0 


01 


20 


01±0 


01 


18 


90±0 


01 


-8 


05±0 


15 


14461^ 


0.75l» 


81 


12 


27 


52 


92 


44 


04 


53 


42 





11 


(2 


0) 





14 


21 


08±0 


01 


20 


42±0 


01 


19 


44±0 


01 


-7 


63±0 


15 


1418li 


0.59lg 



1.12 
1.2(1 



+0.98 

1.29 



-10.63±g:i« 



-8.94Z 
-10.53 
-13.34 
-8.03^ 



+0.17 



+0.18 
0.22 



-8.57 



-0.34 
-0.17 



-10.961°;^ 
.14.10+0-24 

-14.881^1? 



-15.301°;" 
-11.491°;^ 
-11.21±8;|1 



a The * indicates clusters in common with 



Gelatt et al 



12001). For these clusters we provide also the ID assigned by 



Gelatt et al 



1 20011 ) in Table 3. The f 



indicates clusters with B — V>0, V — I < (see Section l4.2l l. The * indicates clusters associated with Ha emission as shown in Fig. [3] 
b Intrinsic effective radius in arcsec and parsec (adopting a distance of 3.82 Mpc from A08) derived with ISHAPE (see Scction l3.2l l. 
c Clustcr cllipticity computed as 1 — r), where r\ is the ISHAPE minor over major axis ratio (see Section 13.21 1. 
d Cluster total magnitudes in the ACS Vegamag system (see Section l3.3l l. 

c Absolute Johnson-Cousins V cluster magnitude computed adopting a distance of 3.82 Mpc from AOS and correcting for a Galactic reddening of E(B — V) = 0.019 
(see Section l3.3ll . A reddening correction of E(B — V) = 0.3 was applied to cluster 51 (see Section |4. 21 . 



Cluster age and mass derived with the Z=0.004 Padova models IGirardi et al and assuming a Salpeters' IMF down to 0.1 Mq. Clusters with t were fitted 

using the Z=0.004 GALEV models (see Section |4,2| |, Only Galactic reddening was assumed, with the exception of cluster 51, for which we solved for reddening 
and derived E(B - V) ~ 0.3. 

g Cluster absolute V magnitude at 10 Myr obtained with the Z=0.004 Padova models (see Section |4~4) | . and with the Z=0.004 GALEV models for f . 
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3.2. Sizes and ellipticities 

Intrins ic sizes of th e candidate star clusters were derived using the ISHAPE task in 



BAOLAB flLarsen 



19991 ). ISHAPE models a source as an analytical function convolved 
with the PSF of the image. For each object, ISHAPE starts from an initial value for the 
FWHM, ellipticity, amplitude, and object position. These parameters are then adjusted in 
a x 2 iterative minimization, until the best fit between the observed profiles and the model 
convolved with the PSF is obtained. We adopted a fitting radius of « 6 x FWHMpg F 
in the ISHAPE procedure, corresponding to ~ 0.5" (~ 14 pixels). We implemented 
both a King (1962) profile with a concentration parameter c = 30 (where c is the ratio 
between the tidal radius and the core radius), and a Moffat profile with a power index of 
1.5. Both profiles provide reasona 



McLaughlin fc van der Marel 



2005 



:>le descriptions of the brightness profiles of GCs (see 
for an extensive discussion of GC brightness profiles). 
The ISHAPE output includes the intrinsic major axis FWHM, minor/major axis ratio, x 2 ! > 
flux, and signal-to-noise ratio for each object, plus a residual image. Intrinsic effective radii 
R e were obtained multiplying the intrinsic FWHM, averaged over the major and minor 
axes, times the conversion factors (1.48 for KING30 and 1.13 for MOFFAT15) provided 
in Table 3 of the ISHAPE tutorial. King and Moffat profiles do not provide significant 
difference in the minimum x 2 or i n t ne final intrinsic R e . We provide in Col. 4 of Table 1 
the R e averaged over the F435W, F555W, and F814W images, both in arcsec and in parsec 
adopting a distance of 3.82 Mpc from A08. For some clusters a best fit could not be found 
in one of the three bands. In this cases, the final R e was obtained by averaging only in two 
bands. 

Ellipticities were derived as e = 1 — 77, where 77 is the ISHAPE minor over major axis 
ratio averaged over the three bands. The values for the individual clusters are provided 
in Col. 5 of Table 1. The R e and e distributions for blue (m F555W — m F814 ^ < 0.9 or 
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m F4 35W - m F555W < 0.5) and red (m F555W - m F8 uw > 0.9 and m Fi ^ w - m F555W > 0.5) 
clusters (see Section [3T31 and Table 1) are shown in Fig. R e varies from ~ 1 pc to ~ 9 
pc, with a median value of ~ 3.2 pc, while the ellipticity ranges from e ~ to 0.4, with a 
median value of ~ 0.1, and with the most elliptical clusters (e > 0.3) being all blue. 



3.3. Photometry 



Aperture photometry with the PHOT package within IRAF was performed at the 
cluster positions within a radius aperture R p hot = ^/ FW H Mm 15 + FW HM PSF) where 
FWHM M15 is the ISHAPE intrinsic cluster size obtained with a MOFFAT 15 profile, 
averaged over the major and minor axes and over the three photometric bands, and 
FWHMps F ~ 2.6 pixels (~0.1"). The sky was computed in a 5 pixel annulus at 5xR p hot- 
Aperture corrections from Rphot t° "infinity" were computed from the results of the artificial 
cluster experiments described in Section 13.44 and applied to the photometry to obtain total 
magnitudes. The artificial cluster tests also provide an estimate of the photometric error as 
a function of magnitude, cluster size, and photometric band. Finally, the p hotometry was 



calibr ated into the HST Vegamag system adding the zeropoints provided in 



Sirianni et al. 



fl2005h . The total F435W, F555W, and F814W magnitudes are given in Cols. 6, 7, and 8 



of Table 1. In Col. 9 we pro vide the abso 



with the transf ormations of ISirianni et al. 



E(B-V)=0.019 ( jSchlegel et al. 



u te V magnitudes in Johnson- Cousins obtained 



(120051 ) . corrected for a Galactic-extinction of 



19981 ) . and adopting a distance modulus of 27.91±0.15 



(3.82 ± 0.27 Mpc, A08). I n Fig . [6] we compare the observed V magnitudes with those 



measured by 



Gelatt et al. 



(120011 ) for the 37 clusters in common. The agreement is pretty 
good, apart from some discrepancy at the brightest magnitudes. Bright clusters tend to be 
located in crowded regions of high star formation, and their photometry is more uncertain. 



In particular, there may be some differences in the adopted photometry aperture between 
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us and 



Gelatt et al. 



(j200lf ). since their photometry extends "up to the radius where the 



clusters visually merged into the background'' 



3.4. Artificial cluster experiments 

We ran experiments with artificial clusters to evaluate the level of completeness and 
photometric error for different magnitudes, cluster sizes, and crowding conditions. These 
tests also provide aperture corrections to the cluster photometry from the adopted R p hot 
aperture to "infinity". 

First, we generated models of extended sources using the Baolab task MKCMPPSF 



(ILarsen 



19991 ) . This task creates a model source by convolving a user-supplied PSF (in this 
case the ACS PSF created from our images) with an analytic profile (we chose a MOFFAT15 
profile) with variable FWHM. We simulated model sources with intrinsic FWHM of 2, 3, 4, 
5, 7, 9, 11 and 13 pixels. 

Then, using the MKSYNTH task within Baolab, we generated (2000x2000) pixel 2 
zero-background images with the model sources homogeneously distributed in a 10x10 
array. In this way the artificial cluster separation is 200 pixels. In MKSYNTH, each source 
is built up by assigning to each photon a random position with probability determined 
by the PSF and the intrinsic source profile. The number of electrons N e added for 
an object of magnitude mag is N e = MKSYNTH. ZPOINT x io-°- 4xma f , and thus 
MKSYNTH. ZPOINT = t exp x 10 0AxZP , where t exp is the total exposure time and 
ZP is the photometric zero point in the specific photometric band. The procedure was 
repeated for all the possible combinations of size (FWHM=2, 3, 4, 5, 7, 9, 11, 13 pixels) 
and magnitudes (mnsTyega from 17 to 22.5 mag, with a step of 0.5), resulting in 12x8 
= 96 images per filter. Each synthetic image of given (FWHM, my ega ) was added to the 
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real image in one of the three galaxy subregions shown in Figs. [7] using the IMARITH task 
within IRAF, and the data reduction process was repeated in the same way as for the real 
stars (see A08) in a (2000x2000) pixel 2 portion of the original science image. The three 
selected fields account for the different crowding conditions within the NGC 4449 field of 
view, with Field 1, chosen on the galaxy center, being the most crowded one, and Field 3, 
selected on a low-background region far from the galaxy center, being the least crowded one. 
Clusters were selected from the final photometric catalog, calibrated as the real star catalog, 
adopting a sharpness cut of sharpnessoAo > 0.5 and, for the F814W images, a magnitude 
cut of m F81iW < 22. The sharpness/magnitude selected catalog was then cross-correlated 
with the coordinate list given in input to MKSYNTH to generate the synthetic images. The 
ISHAPE task was run on the matched catalog to determine the size of the artificial clusters. 
Aperture photometry with PHOT was then performed within an aperture radius R p hot in 
the same way as for the real clusters (see Section \3.3\) . For each filer, this process was 
repeated 12 (mags) x 8 (FWHMs) x 3 (fields) times, providing 28,800 simulated clusters 
in each band. The results of the artificial cluster experiments are summarized in Fig. [8] 
and EO In Fig. [HI we show the fraction of recovered clusters as a function of the F555W 
magnitude for clusters of different sizes in the three selected fields. The completeness is 
close to 100% at the brightest magnitudes, and decreases with increasing magnitude. At a 
fixed magnitude, the comp leteness is lower fo r mo re extended objects and for more crowded 



fields, as already found by iMora et al. 



(120071 ) and 



Smith et al 



fl2007h . 



We also notice that stochastic effects induced by the IMF sampling, not included in 
our study, can affect the cluster detecti on rate; however, this effec t starts to be significant 



for clusters with masses below 10 4 M n ( Silva- Villa fc Larsen 



20111 ). while the majority of 



our detected clusters have masses above this limit (see Section l4~2l) . 



Fig. [9] shows the input minus output F555W magnitude versus the input magnitude 
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for the recovered artificial clusters in the three fields. In this example, the clusters have 
an intrinsic FWHM=3 pixel (~0.1"). The input magnitude is the total cluster magnitude, 
while the output magnitude results from the photometry within an R p hot radius aperture 
(see Section [3731 for details). Therefore, the average mag(inp)— mag(out) is the aperture 
correction from R p hot to "infinite" to be applied to the real cluster photometry. We notice 
in Fig. IH] that the aperture correction is the same in all the three fields at the brightest 
magnitudes. Then it decreases as magnitude increases, and this effect is larger as the field 
is more crowded (~0.2 mag difference at F555W=22.5 in Field 1). This trend is due to 
the increasing importance at fainter cluster magnitudes of the contamination from other 
sources falling within R p hot- We determined for each real cluster the "most representative 
field" by computing the average surface brightness in a (30-100) pixel annulus from the 
cluster center, and comparing this value with the average surface brightness of the three 
fields. Then, given the field, the band, and the cluster R p h t magnitude, we applied the 
appropriate correction to "infinity" , which accounts also for the effect of blending with 
other sources. The 1 a dispersion level around the average mag(inp)— mag(out) value was 
adopted as photometric error. 

The artificial cluster experiments provide also a test on the reliability of the cluster size 
measurements. We show in Fig. [TD]the average input minus output FWHM for the artificial 
clusters as a function of the F555W magnitude, for different FWHM input values and for 
the three test fields. Fig. [10] reveals a trend toward recovering larger cluster sizes, with the 



differences in F WHM being 



agreement with iMora et al. 



arger at fainter magnitudes and for more extended objects, in 



( 120071 ). The differences tend to increase from Field 3 to Field 

1, i.e. going from low-background regions to high-background, crowded regions. In Field 

2, we derive an average FWHM(inp)— FWHM(out)~ —0.06 pixel (~ 1% difference) for a 
cluster with a magnitude of F555W=20 and a size of FWHM=5 pixel, while the difference 
is as large as 0.3 pixels (~6 %) at F555W=22.5. 
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4. Results 

4.1. Luminosity Functions 

We show in Fig. [11] the observed and completeness-corrected cluster luminosity 
functions (LFs) in F435W, F555W, and F814W. The completeness-corrected LFs were 
obtained "replicating" 1/ f carnp i times each observed cluster of given position in the 
NGC 4449 field of view, size, and (mF435w, ^F555VK, ttifsuw) magnitudes, where f CO mpi 
is the product of the completeness in the three bands for that cluster size and crowding 
conditions (e.g. Field 1, 2 or 3, see Section I3T4"|) . and varies from to 1. 

The histograms in the top row of panels of Fig. [11] refer to the total sample. We are 
almost 100% complete at magnitudes brighter than mp 435W, mp 555 w ~ 20 and mp^uw ~ 19, 
while at fainter magnitudes we are strongly incomplete, and the correction to the LF is very 
large. At magnitudes fainter than ~22, where our completeness drops to zero, we are not 
able to recover the LFs. 

Particular attention must be payed when deriving the slope of the LF, because the fact 
that the magnitude detection limit of the clusters shifts toward brighter magnitudes with 
increasing clusters size (see Fig. [8]) can create "artificial" slopes. At the faintest magnitude 
bins, where only the most compact clusters are detected, we are not able to recover the 
contribution to the LF from the extended clusters. For this reason we restrict our analysis 
to clusters with R e < 4 pc and to magnitudes brighter than ~21, which is the detection 
limit for clusters of size R e ~ 4 pc. For magnitudes fainter than this we are not able to 
properly correct the LFs. In the middle and bottom panels of Figs. [11] the size-selected LFs 
are shown separately for blue (mp 555 w ~ ^fsuw < 0.9 or m Fi35W — mp^w < 0.5) and red 
(^F555iy — m F8uw > 0.9 and rrip^w ~ m Ftt$w > 0.5) clusters, roughly corresponding to 
young clusters and old globular clusters (see Section 14. 2D . We derived linear least squares 
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fits to the histograms in the form log N = am + b. Then, from the slope a, we derived the 
slope a of the luminosity function dN(L) oc L a dL using a = —(2.5a + 1). The values of a 
derived for the blue and red cluster populations are provided in Fig. [TT] and in Cols. 2 and 
5 of Table [2j For the young clusters, the slopes a are flatter than ~ —1.5. 

These slopes should be taken with caution. In fact, it has been shown that in the low 
number regime, dire ct fitting to binned histograms can lead to systematic errors in the 



derived indices (e.g. 



Maiz Apellaniz fc Ubeda 



20051 ) . One way around this problem is to 



use a maximum likelihood estimator (MLE), as shown by iMaschberger fc Kroupal (120091 ). 



Following this approach, and using the same blue and red cluster samples as for the binned 
histogram fitting (BHF) method, we derived new slopes, which we provide in Cols. 3 and 
6 of Table [2j The BHF and MLE approaches provide consistent results, within the errors. 
In Table [2] we also report the magnitude interval adopted to derive the slopes in each band 
and for the blue and red cluster samples, respectively. 

It is interesting to compare our slopes with those derived in the literatu re for other 



Gielesl fl2010h 



galaxies. In Fig. [T2]we present an updated version of the diagram proposed by 
which shows the results of power law fits to LFs of clusters in different galaxies taken from 
the literature. In the V band, the slopes fall in the range —2.7 < a < —1.4. Spiral galaxies 
exhibit slopes typically steeper than a ~ —2. If we exclude blue compact galaxies (BCGs), 



there 


is a trend 


Gieles 


2010 


, for 



, for a discussion). BCGs deviate from this trend exhibiting flat slopes at bright 



luminosities. Indeed it has been noticed that the cluster 



than the typical values found in spirals (lOstlin et al. 



2003 



in BCGs tend to 



3e flatter 



Adamo et al. 



2010 



2011aU bh. 



A possible explanation is that the environment in these starburst galaxies has favored the 
formation of massive clusters and /or the disrupt ion of low mass clusters, even if blending 



can flatten the observed LF (I Adamo et al. 



2011ah . The slopes derived in NGC 4449 tend to 
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be flatter than the values found in spirals, and more similar to those of BCGs; its location 
in the a versus luminosity diagram appears intermediate between that of spirals and BCGs. 



4.2. Ages and Masses 



The F435W-F555W (B-V) versus F555W-F814W (V-I) color-color diagram for the 
candidate star clusters is shown in Fig. [13j We can distinguish four different populations in 
the diagram: red clusters (likely old globular clusters) with V — I > 0.9 and B — V > 0.5; 
blue (and young) clusters with < V — I < 0.9 and < B — V < 0.5; very blue (and very 
young) clusters with V — I < and B — V < 0; V-bright cluster s, with V — I < a nd 



Fedotov et al 



20111). We 



B — V > (likely young clusters embedded in ionized g as, see e.g. 
count 4 very blue clusters (clusters number 13, 55, 57, and 73), 6 V-bright clusters (clusters 
number 10, 12, 49, 51, 56, 78), 44 blue clusters, and 27 red clusters. In the figure, we also 
indicate the 11 clusters of Fig. |3] associated with Ha emission. Nebular emission is present 
in all the 6 V-bright clusters, in 2 of the very blue clusters, and in 3 of the blue clusters. 



We overplot the Padova-20 10 simple stellar population (SS 



2010) and the GALEV models (jAnders fc Fritze-v. Alvensleben 



) m odels^ 



initial mass function (IMF), a Galactic reddening E(B-V)=0.019 (ISchlegel et al. 



2003) for a 



Girardi et al. 



Saloeteri ( 19551) 



met a 



licities from Z=0.0004 to Z=0.02. The Padova SSPs are based on the 



1998) and 



( 120081 ) models (e.g., the 



Girardi et al 



P-A GB modeling from 



(2000) 



Marigo et al 



trac ks up to the early-AG B + the detaile d 



Marigo & Girardi! (120071) for M < 7M Q + the 



( 119941 ) models for M > 7M Q ) with the iGirardi et al.l ( 120101 ) correction for low-mass 



Bertelli et al. 



lo w-met allicity AGB tr acks. The spectr a 



m 



Girardi et al. 



(120081 ) and 



Loidl et al. 



libra ry and bolometric corrections are described 



(1200 ll ). The GALEV models are based on the 



4 downloaded at http / /stev . oapd . inaf . it/cgi-bin/ cmd 
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( I1997L 



i sochrqnes o f iGirardi et al.l (120001 ) and on the model atmosphere spectra from lLejeune et al. 



1998), and include gaseous continuum and line emission from ionized gas. The effect 
of the nebular emission on the integrated colors is visible as the turn-over of the models at 
the youngest ages. 

We notice that the Padova models account for the red, blue, and very blue cluster 
populations, but can not account for the 6 clusters with V — I < 0, B — V > 0. If we 
exclude cluster 56, in all the other 5 clusters the nebular emission is spatially coincident 
with the cluster stars (see Fig. [3]), and thus the integrated colors are highly contaminated 
by the ionized gas. This is why the GALEV models, which include a treatment of the 
nebular emission, provide a better match than P adova's. According to the Ha-morphology 



classification proposed by 



Whitmore et al. 



(120111 ) these are emerging clusters (category 3) 



with a mean age of ~3 Myr. On the other hand, the GALEV models fail in reproducing 
the very blue clusters with V — I < and B — V < 0, for which the Padova models provide 
a good match and indicate ages of a few Myr. Two (13, 55) of these clusters are shown 
in Fig. [3j the Ho; emission is not coincident with the cluster itself, but has the form of an 
envelope surrounding the cluster, likely pushed out by massive star winds and supernovae. 
Thus the cluster integrated colors, measured on a relatively small aperture, are not affected 
by the ionized gas c ontrib ution, and are well reproduced by the Padova models. According 



to 



Whitmore et al. 



(120111 ). these clusters belong to category 4a, with a mean age of ~5 Myr. 



Models of different metallicities are strongly degenerate in the (B— V) versus (V— I) 
diagram. In fact, while we can safely exclude a metallicity as low as Z=0.0004, which 
fails to reproduce the red cluster population, models of higher metallicity (Z=0.004, 0.008, 
and 0.02) are all consistent with the cluster colors. We notice in particular that the 
Z=0.02 models are those that best reproduce the blue cluster population, and this holds 
in particular for the GALEV models. This is surprising, because abundance estimates in 



- 24 - 



(Talent 


1980: 


Hunter et al. 


1982: 


Martin 


1997) 



of 12 + logf O/H) JV4449 ~ 8.31; adopting the solar abundance determinations by 



Caffau et al. 



(12008 



20091 ) of 12 + log(O/H) ^ 8.76 and Z Q w 0.0156, and assuming that the chemical 
composition of NGC 4449 follows the solar partition, this translates into Zn4M9 = 0.0055. 
Furthermore, a metallicity of Z~0.004 is consistent with the colors of the bulk of the red 
giant branch (RGB) in NGC 4449 (see A08). Thus, the models that seem to best reproduce 
the blue cluster populations have a metallicity that is ~ 3 .6 times high e r than that of the 



Gelatt et al. 



(120011 ) noticed that 



HII regions and of the RGB stars in NGC 4449. Similarly, 
the ensemble of clusters in NGC 4449 are fitted better with cluster evolutionary tracks 
of somewhat higher metallicity than one would have predicted from the nebular oxygen 
abundance. 

Other sources of uncertainties are the extinction toward the individual clusters and 
the IMF. The reddening vector shown in the panels of Fig. [13] runs parallel to the age<l 
Gyr cluster sequence (for the GALEV models this does not hold for ages younger than 
~10 Myr), implying that an age-extinction degeneracy is present for the blue clusters. 
However, the bluest clusters with B — V < 0, V — I < 0, which fall in proximity of the 
bluest Padova models, are not compatible with a high internal extinction. For the GALEV 
models, we show the effect of reddening in the left panel of Fig. [TH The blue clusters are 
consistent with ages younger than ~100 Myr if the reddening is sufficiently high (up to 
E(B - V) ~ 0.2). The colors of the red clusters (V — I > 0.9, B — V > 0.5) imply old ages 
(>1 Gyr) and negligible internal extinction. Two of the V — I < 0, B — V > clusters 
are compatible with negligible extinction (49, 56), one implies E(B — V) ~ 0.3 (51), while 
three (10, 12, 78) are too blue in V — I for their (B — V) at any reddening. It is actually 
possible that these are not clusters but background emission-line galaxies. 



The effect of the IMF is shown in the right panel of Fig. [TH for the GALEV models, 
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where w e plotted, fo r a Z=0.004 metallicity, Salpeter' s, Kroupa' s (jKroupall200ll ). and 



Scalo' s (1 Scalo 



19861 ) IMFs. A Salpeter and Kroupa's IMF provide similar integrated colors, 
while some difference is observed at the youngest ages for a Scalo IMF. In general, the effect 
of the IMF on the integrated colors is very small at ages older than ~ 1 Gyr. 

For the three metallicities Z=0.004, 0.008, and 0.02, we derived the cluster ages 
minimizing the difference between models and data colors (B— V, V— I, and B— I) according 
to a x 2 criterion, varying the age between ~4 Myr and ~13 Gyr. The model colors were 
corrected for a Galactic reddening of E(B-V)=0.019. We did not solve for possible internal 

ative ly short wavelength 



extinction, because of the small number of bands (3) and the r e. 

range covered by our data. Following the prescriptions of lAvnil ( 119761 ). the 1 a errors were 
derived as the limits of the 0.68 confidence level region defined by x 2 < Xmin + 1> where 
X 2 is the actual value and not the reduced one. Cluster masses were obtained from the 



I banco 



mass-to-light ratios of the best-fitting models in the I banco The procedure was applied 
using both the GALEV models and the Padova models. When using the Padova models, 
we did not fit the V — I < 0, B — V > clusters. For these clusters, the best fitting solution 
was searched with the GALEV models solving for reddening too, since at V — I < the 
GALEV models allow for a disentanglement of the age-reddening degeneracy. For clusters 
10, 12, 78, with B — V > 1.3, a satisfactory fit could not be obtained. 



The results are shown in Fig. [15] in a mass versus age plot for the three assumed 
metallicities. We have adopted a Salpeter IMF between 0.1 and 120 M Q and a Galactic 
extinction E(B-V)=0.019. The age solutions obtained with the Z=0.004 models span the 
whole Hubble time, and the masses are in the range ~ 10 3 — 2 x 1O 6 M . We notice 
particular concentrations of clusters in correspondence of some age bins, such as at ~ 6 



5 Notice that an IMF flatter than Salpe ter's below ~ 1M Q , as suggested by different 



studies (see e.g. Fig.2 in 



Bastianetal 



2010l ). results in cluster masses ~30%-40% lower. 
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Myr, 250 Myr, and 1.6 Gyr. However, as discussed in iBik et al.l ( I2003af ). these peaks are 
not physical, and result from the difficulty in deriving a good fit to the data in those regions 
were the model colors change very rapidly with time. Larger metallicities tend to provide 
younger ages, because of the age-metallicity degeneracy, and lower cluster masses. The 
Z=0.02 models provide a bimodal age distribution peaked at ~ 6 Myr and ~1.5 Gyr. At all 
metallicities, the GALEV models tend to produce younger ages than the Padova models. 

In Fig. [15] we see a trend in the sense that the cluster mass tends to be higher with 
the cluster age. This is caused by a combination of effects. The lower envelope of the 
distribution is a completeness effect: because of the fading of the stellar populations with 
age, old low-mass clusters are not detected. To better show this effect we plot in Fig. [15] 
the line of a cluster of magnitude mp^w = 21.5 at different ages at the distance of 
NGC 4449. The upper envelope in the log(mass) versus log(age) distribution is instead 
due to a sampling effect. In fact, in a log scale, bins of equal size sample age intervals 
which increase with age. The probability of sampling a high mass cluster decreases with 
decreasing age interval, and this is why we do not observe high mass clusters at the youngest 
ages. See however Section [5] for more details. We derive a total mass in clusters of ~1.7, 
~1.1, and ~0.9 x 10 7 M© assuming metallicities of Z=0.004, 0.008, and 0.02, respectively 
(Padova models). These values are lower limits since low- mass old clusters are lost due to 
incompleteness. 

We provide in Cols. 10 and 11 of Table 1 the cluster ages and masses derived with 
the Z=0.004 Padova models, since this metallicity is the closest to the nebular abundance 
and to the metallicity of the RGB stars in NGC 4449, even though we are not able to 
discriminate between the different metallicities from the integrated cluster colors alone. 
For the V — I < 0, B — V > clusters, the solution, when provided, was obtained with 
the GALEV models. As it can be seen in Table 1, the ages of the clusters with associated 
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Ha emission, as derived from the integrated B V, I colors, tend to fall in the range ~ 
5-7 Myr (clusters 13, 49, 51, 55, 56); clusters 9 and 43 have ages of ~ 100 Myr, while for 
cluster 42 we derive an age as old as ~ 300 Myr. We notice in Fig. [3] that for the youngest 
clusters the Ha emission is closely associated with the cluster itself and has the form of an 
envelope surrounding it, while in the oldest ones the nebular emission is more extended 



and/or distant from the cluster, likely because it was pushed out 
supernovae. These findings are in agreement with the results of 



jy massive star winds a nd 



Whitmore et al. 



(2011) 



pointing toward a close connection between the Ha morphology/extent and the cluster age. 

We analyze how clusters of different ages are distributed within the NGC 4449 field 
of view, and how they compare with stars in different age ranges as selected by A08. In 
Fig. [16] we show the spatial distribution of stars and clusters in four age bins (<10 Myr, 
10—100 Myr, 100—1000 Myr, and >lGyr). For the clusters, the age selection is based on 
the results obtained assuming Z=0.004, for a self-consistent comparison with the selection 
of A08. There is a tendency for the distributions of stars and clusters to follow each-other. 
Young clusters are preferentially located in regions of young star formation. Old clusters 
are distributed over the whole NGC 4449 field of view, like the old stars. The difference 
between the spatial distribution of young clust ers/stars and older s tellar populations has 



also been not iced and quantified in the LMC ( 



galaxies (e.g. 



Annibali et al. 



2003 



Bastian et al. 



Bastian et al 



2009bl ) and in other dwarf 



20111 ). However, it is intriguing that in 



NGC 4449 the old clusters do not appear to be distributed homogeneously. As indicated in 
Fig. [lEl a number of linear structures are evident. It is possible that this distribution is a 
reflection of past satellite accretions in NGC 4449, possibly linked to its current ongoing 
star burst. 



Reines et al. 



(120081 ) identified in NGC 4449 13 thermal radio sources, likely embedded 



massive clusters, and derived their physical properties from the nebular emission of the HII 
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regions and from spectral energy distribution fitting to HST UV-near infrared data. The 
radio-detected clusters have ages <5 Myr and follow the distribution of the stars younger 
than 10 Myr in Fig. [TJ3 Three of their thermal radio sources are in common with our 
cluster list, namely cluster 12, 51, and 55, which are displayed in Fig. [3j For clusters 51 and 
55 we derive ages of 5 ± 1 Myr and 6 ± 1 Myr, respectively. 



4.3. Derivation of T 



Several authors have tried to estimate the total mass of clusters recently formed, and 
to compare this mass with the star formatio n rate oyer th e same time scale. The ratio of 



20081 1 . For ages younger than ~10 



stars forming in clusters is referred to as T (iBastianl 
Myr and masses above ~ 1O 4 M , our cluster sample is likely to be complete (excluding 
young embedded clusters). Thus we computed the total mass in clusters more massive 
than this limit and younger than ~10 Myr, and extrapolated it down to 1O 3 M assuming a 



cluster initia 



Adamo et al. 



mass function with slope a = — 2 (a similar approach was adopted by, e.g., 



2011bl ). to obtain the total mass in recently formed clusters. This mass was 



then divided by 5 Myr to obtain the recent cluster formation rate (CFR); younger clusters 
are likely to be lost in our images because still embedded or highly obscured. Considering 
all the results obtained with the Padova and GALEV models for different metallicities, 
we obtain CFR = 1.2 ± 1.9 x lQ~ 2 M^y r~ l which provides , adopting a cu r rent S FR of 



lM & yr 1 from 



McQuinn et al. 



(1201(1 ). T ~ 0.01 ±0.02. 



Goddard et al. 



(120101 ) found 



a positive correlation between V and the star formation rate density T^sfr for a sample 
including spiral and dwarf irregular galaxies. In Fig. [T7] we present an updated version of 
the Goddard et al. diagram in which we include newer results from the literature, as well as 
the values derived for NGC 4449. Here, YIsfr ~ 0.04M & yr~ 1 kpc~ 2 was obtained dividing 
the current SFR by the total field of view of the two ACS pointings. Fig. [17] shows that 
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NGC 4449 lies below the Goddard relation, exhibiting a too low CFR for its current SFR. 
However, if we include the central SSC, with an age between 6—10 Myr and a lower limit 



for the mass of 4 x 1O 5 M (IBoker et al 



20011 ) . T becomes >0.09, and NGC 4449 is in better 



agreement with the re 
the masses derived by 



ation. Another way to compute the CFR in NGC 4449 is to adopt 



Reinesetal 



( 120081 ) for radio-detected embedded clusters, with ages 
<5 Myr. From the masses listed in their Table 5 and adopting an age range of ~5 Myr, we 
obtain CFR = (4.3 ± 0.5) x 10~ 2 M Q yr- 1 and F ~ 0.04 ± 0.01, higher than the previous 



value without the central SSC, but still below the Goddard relation. 



4.4. Cluster classification 



The terms of "super star cluster ", "populous clust e r", or "massive cluster" are not 
univocally defined in the literature. lLarsen fc Richtlerl ( 120001 ) defined young massive 
clusters (YMCs) to be those with (B -V) < 0.45 and My < -8.5, for (U - B) > -0.4, 
or M v < -9.5, for (U - B) < -0.4. We do not have ACS imaging of NGC 4449 
in the U band, but from the SSP models we see that a (U — B) = —0.4 roughly 
corresponds to m Fi3 5w — m F555\v ~ 0.09, accounting also for the Galactic reddening toward 
NGC 4449. Thus we can select YMCs to be those with (0.09 < m Fi35W — m F555W < 0.45, 
M v < -8.5), and (m F435W - m F555W < 0.09, My < -9.5). We end up with 12 YMCs, 



plus the central SSC. This provides a cluster specific frequency of T, 



iv 



0.65 



M r( N GC 4449) = —18.2. For their sample of 21 nearby spiral galaxies, 



, given 



Larsen fc Richtler 



( 120001 ) derived Tm in the range 0.00—1.77, with an average value of ~0.60 and a median 



Larsen k Richtlerl (I1999[ ) to be 



6 The young cluster specific frequency was defined by 
T N = N x io°- 4x ( M s+ 15 ), where N is the number of YMCs, and M B is the total absolute B 
magnitude of the galaxy. 
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of 0.44. Thus, the specific frequency of YMCs in NGC 4449 is above the median and 
average values derived in spirals, likely because of its high SF activity, but not as extreme 
as the frequency derived in the most rich spirals. For comparison, the LMC has T/v = 0.57 



(ILarsen fc Richtler 



20001 ) . What is the cluster specific fr equency in starburs ting irregulars 



Billett et al 



(20 021) and A08 we 



such as NGC 1705 and NGC 1569? From the results of 

have i n NGC 1705 only 2 clusters (including the SSC) satisfying the ILarsen fc Richtler 
Jl999h definition ; this implies T N ~ 0.87, adopting Mb {NGC 1705) = -15.9. From Table 2 



of 



Hunter et al. 



(120001 ) . and adopting the new distance of 2.96 Mpc by Grocholski et al. 



(in preparation), we have 23 YMCs in NGC 1569, which imply a frequ ency T/y ~ 1.81 



Larsen fc Richtler 



(Mr ( 1569) ~ —17.76), higher than those of the most rich spirals in the 
(120001 ) sample. Hence, NGC 4449 is more active than average in forming young massive 
clusters, but definitely less than other starburst dwarfs. 



Gelatt et al. 



(]200ll ) introduced a new definition of "populous clusters" and "super star 



clusters" to b e tho se with —10 < My < — 11 and My < — 11 at 10 Myr, respectively, while 



Billett et al. 



(120021 ) adopted less stringent limits of —9.5 and —10.5 in My at 10 Myr. These 
definitions are somewhat model dependent, since the age correction to 10 Myr depends on 
the particular model used. We provide in Tab le [3] the number of populous and super star 



clusters identified with the 



Gelatt et al. 



(120011 ) definition, assuming different metallicities 
(Padova models). For each metallicity, the range in numbers reflects the uncertainty in age. 
The My at 10 Myr is given in Col. 12 of Table 1. We identify 4 to 12 young (< 1 Gyr) 
SSCs, and a conspicuous population of old (>1 Gyr) SSCs (21 to 27) that could be massive 
globular clusters. The Z=0.02 metallicity provides the lowest number of SSCs. Assuming a 
metallicity of Z=0.004 or Z=0.008, we end up with a mini mum number of 7 young SSCs. 



This is larger than what was found by 



Gelatt et al. 



fj200lh . who identified in NGC 4449 



three young SSCs in addition to the central one (their cluster IDs 31, 27 and 47). These 
three clusters correspond to our clusters 43, 48, and 18, which we classify as young SSCs 



31 



as well. Cluster 43 is located just ~ 20 pc east from the central SSC, while cluster 48 is ~ 
150 pc west. Cluster 18 is located N-E at ~1.3 kpc from the galaxy center, in the region of 
young SF. Our additional four young SSCs are cluster 6, fairly close to cluster 18; cluster 
36, which exhibits a high ellipticity and is located ~ 250 pc east of the central SSC; cluster 
66 and cluster 67, ~ 900 pc S-W and ~1.1 kpc south from the galaxy center, in proximity 
of the stream of blue stars extending south. Thus all young SSCs are located in regions of 
active star formation. The old SSCs are instead more homogeneously distributed over the 
NGC 4449 field of view, following the same distribution of the old stars in Fig. [161 Several 
populous clusters, from 8 to 22, are also identified in our images. 



4.5. Correlations between cluster properties. 

We investigated the presence of correlations between the cluster properties. We find 
significant correlations only for the red (V — I > 0.9, B — V > 0.5) cluster population 
between the effective radius R e and magnitude, R e and cluster mass, and ellipticity e and 
magnitude (see Fig. [TBI . The effective radius tends to increase with the cluster magnitude, 
meaning that brighter clusters are more compact, while fainter clusters tend to be more 
extended. The Spearman correlation coefficient is r s = 0.63 with 25 degrees of freedom, 
implying a probability P~0.0004 that a correlation is not present. The correlation between 
R e and magnitude seems, at least in part, driven by mass: in fact we find that more massive 
clusters are also more compact. With the masses derived from the Z=0.004 Padova models, 
the correlation coefficient is r s = —0.54, implying a probability P~0.004 that a correlation 
is not present. Our conclusions do not change if we adopt the masses derived with the 
GALEV models, or with different metallicities. Finally, we find an anti-correlation between 
magnitude and ellipticity, in the sense that brighter clusters tend to be more elliptical. The 
correlation coefficient is r s = —0.49, implying a probability P~0.009 that a correlation is 
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not present. 



In an early study, 



van den Berghl (119831 ) noticed that in the LMC more luminous 



clusters tend to be more flattened, but the data now available show no evidence for 



such a correlation (Ivan den Bergh 



20081 ) . This result is true both for g l obula r clusters 



van den Berghl (120081 ) pointed 



and for younger clusters. For what concerns the SMC, 
out that available data are too scant to determine if there is a correlation between the 
luminosity and the flattening of clusters, but noticed that the four brightest SMC clusters 
are all very flattened having < e >= 0.26. Interestingly, wCen, the brightest globular 
cluster of the MW, is also one of its most elongated ones (see the Harris catalog at 



http : //www. physics .mcmaster . ca/~harris/mwgc . dat ). 



5. Cluster disruption and infant mortality 



The term "infant mortality" was initially introduced by 



Lada fc Ladal (120031 ) to 



indicate that 90% of the embedded clusters in the solar neighbourhood do not survive the 



gas expulsion ph ase. According to theory, th e effects of g as ex pulsi on should be large 



within ~20 Myr (iGoodwin &: Bastian 



20061 ). 



Fall et al 



( I2005h and 



Whitmore et al 



y ove r 



(120071) 



introduced the concept of "long-term infant mortality" , finding that in the Antennae 90% 
of the clusters, independent of their mass, are disrupted every dex of age, a process that 
continues between a few Myr to ~ 100 Myr, and possibly up to ~1 Gyr. The same fraction 
(90%) of infant morta lity was clai med to be present also in other g alaxies, such as the SM C 



Chandar et al. 



Pellerin et al 



2006), the LMC (IChandar et alJl2010a[ ). and M83 (IChandar et al.ll2010ri ). 



( 120071 ) claimed significant infant mortality in the spiral NGC 1313 from 



the fact that O stars were more "clustered" together, over ~100 pc, than B stars (notice 
however that these scales imply that this is unlikely to be caused by traditional "infant 
mortality"). These results are in contradiction with other studies, which instead do not find 
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evidence for long-term infant mortality (IGieles fc Bastianll2008l ; iBastian et al.ll2009al ) 



The most common way to study cluster disruption is to u se the cluster age d istrib ution 



which, however, can be heavily affected by incompleteness. 



Gieles & Bastianl (120081 ) 



proposed a new method of studying cluster disruption b ased on the relation betw een the 



most massive cluster, M max , and the age range sampled. 



that, assuming that clus ters are stochastica 



Gieles fc Bastianl (120081 ) showed 



function with index —2 (jZhang &: Fall 



1999 



V samplec 


from 


Bik et al. 


2003b 



power-law cluster in itial mass 



de Grijs et al. 



20031 ). and that 



the cluster formation rate (CFR) is constant, M max scales with the age range sampled, such 
that the slope in a \og(M max ) vs. log(age) plot is equal to unity. In the scenario of mass 
independent (MID) cluster disruption, the slope turns out to be [1 + log(l — /mid)], where 
Imid is the fraction of clusters destroyed in each age dex. Thus the slope becomes flatter 
than 1 as cluster disruption is more important. 

We investigated the \og(M max ) vs. log(age) relation for the young (age<l Gyr) clusters 
in NGC 4449. Our results are shown in Fig. [19j We adopted an age bin of 1 dex starting 



from log(age[yr])=6. This is the same bin adopted 



than the 0.5 dex bin adopted by 



by. 



Chandar et al 



(l2010bf ). and is larger 



Gieles fc Bastianl (120081 ) . Given the bin size, we will not be 



able to characterize the "short term infant mortality" within the first ~10 Myr. The slope 
of the \og(M max ) vs. log(age) relation was computed for the three metallicities Z=0.004, 
0.008 and 0.02, with the GALEV and with the Padova models. We assumed both no 
internal reddening (top panels), and an E(B — V) as high as .2 (middle panels). Notice 



Hunter fc Gallagherl (Il997f ) 



that this value is twice the average E(B — V) ~ 0.1 derived by 
in NGC 4449 HII regions, and should thus be considered an upper limit. The slopes exhibit 
some dependence on metallicity and on the adopted models (Padova or GALEV). If we 
average over the three metallicities and over the two models, we obtain 0.76 ± 0.07 in the 
case of no internal reddening, and 0.58 ± 0.12 in the case of E(B — V) — 0.2. These slopes 
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imply cluster disruption fractions of Jmid = 0.421qjo and /mid = 0.62i j2, significantly 
lower than the 90% long term infant mortality found by other authors. 

Our conclusions rest on the assumption that the CFR is constant. However, what if 
this is not the case? We can ass ume that the CFR follows the SFR in the field. From Fig. 9 



and 10 of 



McQuinn et al. 



(120101 ) we derive SFRs of ~0.98, -0.36, and ~0.19 M Q /yr in the 
<10 Myr, (10-100) Myr, and (100-1000) Myr age bins, meaning that NGC 4449 has been 
more active at recent epochs than in the past. We thus "correct" the M max values for the 
fact that the SFR was not constant over the past 1 Gyr. Since M mnx o c Nhi. n , where Nu n is 



Gieles fc Bastianl ( 120081 )). and N bin oc 



the number of clusters in each age bin (see Eq.(7) of 
CFR oc SFR , we re-compute the slopes using the normalized log(M max /CsFR), where Csfr 
is proportional to the SFR in the corresponding age bin (bottom panels of Fig. IT9"]) . In 
this case the average slope is as steep as 1.12 ± 0.07, as a consequence of the more intense 
average SF experienced by NGC 4449 over the last ~100 Myr, and in particular over the 
last ~10 Myr. This is consistent with no cluster disruption. 

Finally, we should c omment on the fact that our catalog does not include the central 



ssc. 



Gelatt et al. 



(200 1J) deriyed for the SSC an age between 8—15 Myr and a mass of 



2 x 1O 5 M , while boker et al. 



(120011 ) derived an age between 6—10 Myr and a lower limit 
for the mass of 4 x 10 5 M Q . Thus, if i ncluded in our M m „. x analysis, the SSC would provide 



significantly flatter slopes. However, 



Gelatt et al. 



(l200lh showed the the central SSC in 



NGC 4449 is a very peculiar object: it consists of a bright condensation 0.1" in diameter 
sitting in the middle of an elongated, bar-like structure. The minor-to-major axis ratio o 
the st ructure is 0.6, while the cluster is essentially round in the outer parts. 



Gelatt et al. 



( 120011 ) suggested that this structure is the result of the interaction scenario that also 



produced the counterrotating gas systems and extended gas filaments and streamers in 
NGC 4449. In particular, the peculiar structure at the center of the SSC could be the debris 
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from cannibalism of a smaller companion that has fallen into the center of the galaxy. In 
this context, the central SSC in NGC 4449 appears a peculiar object whose formation is not 
comparable to the formation of the population of "normal" clusters. More generally, it has 



been well documented that nuclear star clusters i n galaxies ( "nuclei" 



van der Marel et al. 



tend to have distinct 



2007, and references 



properties from globular clusters in general (e.g. 
therein). This is not unexpected, given the special location at the center of a gravitational 
potential well. This further motivates our decision to exclude the central SSC in NGC 4449 
from our analysis of the general cluster population. 



Results and discussion 



We have presented a study of 81 candidate star c lusters in NGC 4449, 37 of which in 



common with a previous study by iGelatt et al 



(1200 lh . based on our ACS/WFC F435W 



(~B), F555W (~V), F814W (~I), and F658N (Ha) images. Our ACS images are saturated 
in the center of the galaxy, preventi ng an analysis of the central super st ar cluster (SSC), 



already studied by others authors (IGelatt et al. 



2001 



Boker et al. 



200l|). The spatial 



distribution of the clusters appears to follow that of the resolved stars in NGC 4449. 

We derived the cluster intrinsic properties such as the intrinsic effective radius R e and 
the ellipticity e, as well as the integrated magnitudes. R e is found to vary from ~ 1 pc to 
~ 9 pc, with a median value of ~ 3.2 pc, while the ellipticity ranges from e ~ to 0.4, 
with a median value of ~ 0.1, and with the most elliptical clusters (e > 0.3) being all blue 
(B — V < 0.5 or V — I < 0.9). Artificial cluster experiments show that the errors on the size 
measurements increase from bright to faint and from compact to more extended clusters. 
The error in size is ~1% for a cluster with R e ~4 pc and V=20, and it is as high as ~ 6% 
if V — 22.5. These should be considered lower limits, since they were derived fitting the 
artificial cluster profiles with the same law assumed to create them. Very extended clusters 



-36 - 



(R e ~ 8 pc) suffer strong incompleteness at relatively bright magnitudes, and are easily lost 
in the most crowded regions. 

In a B — V vs. V — I color-color diagram, the candidate star clusters define four 
different populations, a red one with B — V > 0.5 and V — I > 0.9 (27 clusters), a blue 
one with < B — V < 0.5 and < V — I < 0.9 (44 clusters), a very blue one with 
B — V < 0, V — / < (4 clusters), and a population of V-bright clusters with V — I < 
and B — V > (6 clusters). We find nebular emission associated with 11 clusters out of 81, 
and specifically with all the 6 V-bright clusters, with 2 of the very blue clusters, and with 
3 of the blue clusters. SSP models of different metallicities are strongly degenerate in the 
color-color diagram, thus our results are strongly affected by the age-metallicity degeneracy. 
However, independent of metallicity, red clusters have ages > 1 Gyr, blue cluster < 1 Gyr, 
and very blue clusters are as young as a few Myr. Likely, the V-bright objects are young 
(age< 5 Myr) massive clusters embedded in the ionized gas, as suggested by the fact that 
the Ha emission is coincident with the cluster itself, even if we can not exclude that some 
of them are background emission-line galaxies. From the B, V, I integrated colors, we are 
able to derive an age for 8 clusters out of the 11 associated with nebular emission. The 
majority of them have ages of ~ 5-7 Myr, two have ages of ~ 100 Myr, and one is as old 
as ~300 Myr. For the youngest clusters, the Ha emission is closely associated with the 
cluster itself and has the form of an envelope surrounding it, while in the oldest ones the 
nebular emission is more extended and/or distant from the cluster, likely because it was 
pushed out b y massive star w i nds an d supernovae. These findings are in agreement with 



the results of 



Whitmore et al. 



(120111 ) pointing toward a close connection between the Ha 



morphology/extent and the cluster age. 



Notice that clusters as young as a few Myr are numerous i n NGC 15 69 (Hu nter et al 



2000 



Origlia et al 



20011 ) . and are possibly present in IZw 18 (jContreras al 



20111 ), but 
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are not found in NGC 1705 (IBillett et al.l l2002t lAnnibali et all 120091 ) . If we assume 
a Z=0.004 metallicity (which is consistent with the HII region abundances and with 
the RGB colors of the resolved stars in NGC 4449), we find that the cluster ages are 
distributed quite continuously over a Hubble time, and the clusters have masses in the 
range ~ 10 3 — 2 x 1O 6 M . With the the Z=0.02 models we obtain instead a completely 
different scenario, i.e. a bimodal age distribution strongly peaked at ~ 6 Myr and ~ 1.5 
Gyr. This metallicity is a factor ~ 4 highe r than the me t allicit y derived in the HII regions, 



and thus difficult to accept. However, also 



Gelatt et al. 



(1200 ll ) noticed that the ensemble 



of clusters in NGC 4449 are fitted better with cluster evolutionary tracks of somewhat 
higher metallicity than one would have predicted from the nebular oxygen abundance. If 
we consider the results obtained with the Z=0.004 SSP models, we find that, within the 
uncertainties, the cluster age distribution is consistent with the SFH of NGC 4449 derived 
by McQuinn et al (2010) from the analysis of the CMDs of the galaxy resolved stars. 
Furthermore, young clusters are preferentially located in regions of young star formation, 
while old clusters are distributed over the whole NGC 4449 field of view, like the old stars 
(although we notice that some old clusters follow linear structures, possibly a reflection of 
past satellite accretion). 

These results may sound obvious, if one believes that all stars originally form within 
clusters, but actually this i s not always the case . For instance, the well known bimodal 



SFH of the LMC clusters (IPagel k, Tautvaisiene 



of representative LMC regions (see e.g. 



1998 1) is complete ly at odds with that 



Smecker-Hane et al. 



2002 . and Fig. 7 in Tosi 



et al. 2006). Thus, also the scenario of a bimodal cluster age distribution can not be 
excluded a priori. We derive a lower limit of ~ 10 7 M^ for th e total mass in clusters. For 



comparison, the total mass from the SFH by 
view is 3 ± 0.5 x 1O 9 M . 



McQuinn et al. 



3) i 



(2010) m the same field of 
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We investigated possible correlations between the cluster properties (magnitude, size, 
ellipticity, mass and age). Significant correlations are found only for the old (V — I > 0.9, 
B — V > 0.5) cluster population between the cluster effective radius and the magnitude, 
the effective radius and the mass, and between the ellipticity e and the magnitude. Brighter 
and more massi ve clusters tend to be more compact, and brighter clusters tend to be also 



more elliptical. Ivan den Berghl ( 120081 ) noticed that, while in the LMC there is no evidence 
for a correlation between ellipticity and luminosity, the four brightest SMC clusters are all 
very flattened having < e >= 0.26. Interestingly, also cuCen, the brightest globular cluster 
of the MW, is one of its most elongated ones. 

Among the old clusters in our sample, one of the most massive (M ~ 1.7 x 1O 6 M ) and 
elliptical (e ~ 0.2) one (namely, cluster 77), appears surrounded by a symmetric structure 
of blue stars, which could be tidal tails associated with a dwarf galaxy currently disrupted 
by NGC 4449. This cluster is treated in more details in an accompanying paper (Annibali 
et al. in preparation). 

Because of the combined effects of incompleteness, cluster mass-function sampling, 
and possibly mass- dependent cluster mortality (more massive clusters are more difficult 
to destroy), clusters are not distributed uniformly in a log mass versus log age plane, but 
a trend is observed in the sense th at the cluster mas s seem s to increase with the age. 



Following the approach proposed by 



Gieles & Bastianl (120081 ) , we used the upper envelope of 



this distribution, i.e. the maximum cluster mass in different age bins, to test the hypothesis 



of mass-independent cluster disruption. We do n ot find evidence 



long- t erm in fant morality found by other authors (IWhitmore et al. 



;or the high (90%) 



2007 



Chandar et al. 



2006 



2010a 



bj). If we assume that the cluster formation history follows that of the field 



stars our data are consi stent with no cluster disru ption within ~1 Gyr, in agreement with 



Gieles fc Bastianl (J2008I) and 



Bastian et al. 



fl2009ar ). 
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As mentioned in the introduction, NGC 4449 is a very active star forming galaxy, the 
on e with highest star form ation rate in the sample of 18 starburst dwarfs recently examined 



by 



McQuinn et al. 



(120101 ) . T his is likely the cau s e of i ts large number of detectable YMCs, 



confirming the suggestion by 



Larsen fc Richtlerl (120001 ) of a direct correlation between the 



number and luminosity of YMCs and the SF activity of the parent galaxy. We find that 
the YMC specific f reque ncy in NGC 4449 is higher than the average frequency derived by 



Larsen &: Richtlerl (120001 ) in nearby spirals, and also higher than in the LMC, but lower 
than the frequency derived in other starburst dwarfs such as NGC 1705 and NGC 1569. 
These results are confirmed by NGC 4449 young cluster LF (dN(Ly) oc Ly l ^dL) flatter 
than those derived in late-type spirals, and more similar to those of blue compact galaxies, 
perhaps because the environment in these starburst galaxies has favored the formation of 
massive clusters. Furthermore, NGC 4449 seems to host a significant number of SSCs: 
adopting the defini tion of SSC to be a cluster with M v < —11 at 10 Myr, introduced by 



Gelatt et al. 



fj200lh . we find that there are 7 young (age<l Gyr) S SCs in NGC 4449, b esides 



Gelatt et al. 



(j2001[ ) from 



the central one. Three of these clusters were already identified by 
WFPC2 data. All the young SSCs are located near the galaxy center or in regions of active 
SF. 

Of particular interest is the comparison of the cluster properties with the SFR density, 
which is a usef ul physical quan t ity to rank the SF activity of galaxies of different mass 



and size. From 



McQuinn et al. 



( 12010( 1. the average SFR over the last ~100 Myr is ~0.42 
M & /yr, and from Fig. [TBI of this paper we estimate that stars younger than ~100 Myr are 
distributed over an area of ~7.85 kpc 2 . This provides a SFR density of ~0.05 M@/yr/kpc 2 . 
On the other hand, we count 15 clusters (including the central SSC) with ages <100 Myr, 
which imply a cluster density of ~1.9/kpc 2 over the same area. 



For NGC 1705, 



Annibali et al 



(120091 ) give 9 clusters (including the central SSC) 
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younger than ~100 Myr concentrated in the central 26" x 29" area, corresponding, at a 



20011 ) . to ^642x717 pc 2 . This provides a cluster density 



distance of 5.1 Mpc (ITosi et al. 
of 19/kpc 2 . To compute t he SF R density, we sum up the rates in Regions 7, 6 and 5 in 



Table 3 of 



Annibali et al. 



(120091 ) and divide the result by the total area, which gives ~0.08 
M Q /yr/kpc 2 . Thus, despite the fact that the SFR density in NGC 1705 is a factor ~2 
higher than that in NGC 4449, the cluster specific frequency is ~ a factor 10 higher. 

NGC 1569 i s the dwarf galaxy with the highest SFR density: from HST/NIC2 data, 



Angeretti et al 



(120051 ) derived an average SFR of ~ 2.5 M Q yr 1 kpc 2 over the past 100 



normal clusters ( Hunter et al. 



200ol ; 



dc. NGC 156< 


} is ah 


Orielia et al. 


2001) 



so known to host 4 SSCs, and several 



Hunter et al. 



(I2000h identified 48 



clusters in their WFPC2 images, 28 of which fall within the 19" x 19" NIC2 field of view. 
From the cluster photometry in their Table 2, and using the Padova models, we estimate 
that 25 of them have ages <100 Myr. This implies a cluster density of ~357/kpc 2 , and a 
cluster/SFR density ratio of ~142, a factor ~4 higher than in NGC 4449. 



These results agree with the findings of 



Larsen fc Richtlerl (120001 ). according to which 



NGC 1705 and NGC 1569 deviate from the cluster specific luminosity T L (U)-SFR density 
T^sfr relation defined by star forming galaxies, exhibiting very large T L {U) for their T^sfr- 



Another interesting correlation is the one found by 



Larsen! (120021 ) between the brightest 



cluster V-band absolute magnitude, My t9htest , and the galaxy SFR, suggested to be 
caused by size-of-sample effects. This trend has been used to argue for a uni versality 



of clu ster formation, i.e. stochastic sampling from a universal mass function. 



Bastian 



( 120081 ) showed that the absolute magnitude of the brightest young cluster in a galaxy 
is a good tracer of its current SFR; he also showed that the fact that dwarf /irregular 
starburst galaxies often lie signicantly above the relation can be explained in terms 
of their SFH proceeding in short bursts. From our Table 1, the brightest cluster in 
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NGC 4449 has M v = -10.75 ± 0.15, which falls sig nificantly (> 5a) belo w the relation 

M brighte S t = L87 ( ±0 . 6) X l og SFR - 12. 14 



(±0,07) (jWeidner et al. 



2004) . adopting a 



20101 ) . On the other hand, the cen t ral SS C 



Gelatt et al. 



(120011 ) and 



current SFR of ~ lM & yr~ l (iMcQuinn et al. 

with a magnitude of My = —12.54 ± 0.15 (adopting the value from 
correcting for a distance modulus of DM = 27.91 ± 0.15) turns out consistent with the 
relation within 3a. 
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Table 2. Luminosity function slopes 



Blue clusters a Red clusters 

Band a (BHF) c a (MLE) d Fit interval a (BHF) c a (MLE) d Fit interval 



F435W -1.37 ±0.30 -1.36 -11.0 < M Fi3BW < -7.0 -1.60 ±0.35 -1.37 -9.5 < M F435W < -7.0 

F555W -1.47 ±0.22 -1.45 -11.0 < M F555W < -7.0 -1.30 ±0.20 -1.15 -10.0 < M F555M / < -7.5 

F814W -1.47 ±0.12 -1.52 -11.5 < M F81iw < -7-5 -1.10 ±0.22 -1.07 -11.0 < M F814W < -8.5 

a mF555W — rn F8 i4w < 0.9 or mi? 4 3 5l y — m F555W < 0.5. 
b "iF555lv - ™-FSi4W > 0.9 and m F435W — m F555W > 0.5. 
c Binncd histogram fitting. 
d Maximum likelihood estimator. 



Table 3. Cluster classification 



Zssp N cl * JV_io..-n b AU0..-11 b ND_ 10 .._ n c W<-n d 7V<_ n d JVD<_n e 
(<lGyr) (>lGyr) (<1 Gyr) (>1 Gyr) 

[fcpc~ 2 ] [fcpc~ 2 ] 



0.004 81 18-18 1-4 0.7-0.8 7-12 24-27 1.2-1.5 

0.008 81 13-15 3-6 0.6-0.8 7-11 21-24 1.1-1.3 

0.02 81 6-7 2-3 0.3-0.4 4-5 24-25 1.1-1.1 

a Total number of candidate star clusters, not including the central super star cluster which is 
saturated in our ACS images. 

b Number of populous star clusters, taken to be those with — 10 < My < —11 at 10 Myr. Notice 
that the definition is slightly different than that adopted by Billett et al.(2002), Table 2. The range 
in numbers reflects the uncertainty in age. 

c Number density of populous clusters found within the ~ 26 kpc 2 area covered by our ACS data. 

d Number of super-star clusters, taken to be those with My at 10 Myr brighter than —11. Notice 
that the definition is slightly different than that adopted by Billett et al.(2002), Table 2. The range 
in numbers reflects the uncertainty in age. 

e Number density of super-star clusters found within the ~ 26 kpc 2 area covered by our ACS 
data. 
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Fig. 1. — Mosaicked ACS F555W image (380 x 200 arcsec 2 ) with superimpose d the 81 can- 
didate star clusters (red circles). The footprint of the WFPC2 images used by lGelatt et al. 
(120011 ) and their candidate clusters (blue crosses) are shown for comparison. 




Fig. 2. — (4.5x4.5) arcsec 2 F555W images of the cluster candidates. The big cluster right 
of cluster number 43 is the central SSC, which is saturated in our images. 




Fig. 2. — continued 




Fig. 2. — continued 




Fig. 3. — Color-composite images (F435W=blue, F658N=green, F814W=red) of the clusters 
with associated Ha emission. For each cluster, the field of view is (4.5x4.5) arcsec 2 . 
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Fig. 4. — (4.5x4.5) arcsec 2 F814W images of two clusters classified by 
but considered by us to be "stellar associations" (see Section 13. ip . 
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Fig. 5. — Distributions of the intrinsic effective radii R e (top panel) and of the ellipticities 
e (bottom panel) for blue {mFb^bW — mFsnw < 0.9 or mp^w — nT>F555W < 0.5) and red 
("1F555W — "mpsuw > 0.9 and m F435W — m F555W > 0.5) clusters. 
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Fi g. 6. — Comparison of the observed Johnson-Cousins V magnitudes measured by us and 
by iGelatt et al.l (120011 ) for the clusters in common. The solid line is the one-to-one relation. 
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Fig. 7. — Artificial clusters added to the (380 x 200) arcsec 2 mosaicked ACS F555W image 
in three fields of different crowding. From left to right are Field 2 (medium crowding) , Field 
1 (high crowding), and Field 3 (low crowding). 



- 61 - 



l 

0.8 
0.6 
0.4 
0.2 

1 

0.8 

I 0.6 

o 

u 

- 1 0.4 
0.2 

1 

0.8 
0.6 
0.4 
0.2 




FIELD 


1 


FWHM 


= 3 


FWHM 


=4 


FWHM 


= 5 


FWHM 


= 7 


FWHM 


= 9 


FWHM 


= 11 



FIELD 2 


FWHM 


= 3 


FWHM 


=4 


FWHM 


= 5 


FWHM 


= 7 


FWHM 


= 9 


FWHM 


= 11 



\ 



\ 



\ 



H h 



FIELD 3 

FWHM =3 

FWHM=4 

FWHM =5 

FWHM =7 

FWHM =9 

FWHM=11 
I i_ 



\ 

\ 

\ 



_L 



_l L. 



18 



20 

F555W 



22 



inp 



Fig. 8. — Completeness in the F555W band for artificial clusters with a MOFFAT1 5 intrinsic 
profile and FWHM of 3, 4, 5, 7, 9 and 11 pixel (1 pixel=0.035") in the three test fields (the 
crowding increases from Field 3 to Field 1). 
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Fig. 9. — Example of input minus output magnitudes in F555W for artificial clusters with 
an intrinsic M0FFAT15 profile and an intrinsic FWHM=3 pixels (~0.1"). The output 
magnitudes were derived through photometry within a R p hot aperture (see Section 13.31 for 
details), while F555Wi np is the total cluster magnitude. The average mag(inp)— mag(out) 
value at different magnitudes (central red line) is the aperture correction from R p hot to 
"infinite" to be applied to the photometry. The ±1<t levels (blue external lines) quantify the 
photometric errors. 
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Fig. 10. — Average input minus output FWHM (in pixels) for the artificial star clusters as 
a function of the F555W input cluster magnitude, for different FWHM input values and for 
the three test fields. 
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Fig. 11. — Distributions of cluster magnitudes in the three bands (F435W, F555W, F814W). 
The solid histogram is the observed distribution, while the dotted one is the completeness- 
corrected distribution. The histograms in the top row of panels refer to the total cluster 
sample, while in the middle and bottom panles we applied a cut in size of R e <4 pc. The 
separate distributions for blue (rriF555w — rnFsuw < 0.9 or rriF^w — mFbhsw < 0.5) and 
red {mFbbbW — mFsuw > 0.9 and tjif^w — rriF555W > 0.5) clusters are shown. The straight 
lines are the linear least squares fits in the form log iV = am + b to the different distributions, 
obtained considering only the magnitude bins brighter than 21. At magnitudes fainter than 
this we are not able to properly correct the LFs (see Section H~T1 for details). The slope a of 
the fit, and the slope a of the LF in the form dN(L) oc L a dL, are given in each panel. 
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Fig. 12. — Literature luminosity function slopes a in the V band as a function of the 
fit interval for young clusters in different galaxies: spiral galaxies (NGC 62 8, NGC 1313 



NGC 3184, NGC 4395, NGC 5236, NGC 6744, N GC 6946, NGC 7 793 



Mora et al 



Hwang &: Lee 



200flh: the LMC f lLarsen 



20081 ); the Ant ennae (Whit more et al 



ESQ 185-IG13; NGC 1705, 



2002 "); M 51 



Adamo et al. 



2011a. 



Gieles et al. 



2006 



Larsen 



Haas et al 



2002 



2008 



199 91); blue compact galaxies (Mrk930; 



Annibali et al 



20091 ); NGC 4449 (this 



work). For NGC 1705, the leftmost BCG in the diagram, the slope is the average of the 
binned histogram fitting result (a ~ —1.6) and of the maximu m likelihood est i mator result 



la 



—2.3) obtained from the cluster magnitudes reported in 



Annibali et al. 



(120091 ). The 



significant difference between the two values is likely due to the small statistics. 
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Fig. 13.— F435W-F555W (B-V) versus F555W-F814W (V-I) color-color diagram for 
the candidate star clusters. Triangles indicate clusters with asso ciated Ho; emission shown in 
Fig. El Overplott ed are the Padova simple stellar popu 



GALEV models (jAnders fc Fritze-v. AlvensleberJ (120031 ). right panel) for a Salpeter IMF, a 



at ions (IGirardi et al. 



2010h and the 



Galactic reddening E(B-V)=0.019, and different metallicities. For the Z= 0.004 metallicity, 
we indicate the models of log(age[yr])=6.6, 7, 8, 9, and 10. The arrow indicates the reddening 
vector for E(B - V) = 0.2. 
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Fig. 14.— F435W-F555W (B-V) versus F555W-F814W (V-I) color-color diagram for 
the candidate star clusters. Triangles indicate cl usters with associated Ha emiss i on sh own 
in Fig. |3j Overplotted are the GALEV models ( lAnders fc Fritze-v. Alvenslebenl (120031 ) for 
Z=0.004 and for different reddenings (left panel) and IMFs (right panel). Notice that the 
line corresponding to Kroupa' s IMF completely overlaps Salpeter' s line in the right panel. 
For the Z= 0.004 metallicity, we indicate the models of log(age[yr])=6.6, 7, 8, 9, and 10. 
The arrow indicates the reddening vector for E(B — V) = 0.2. 
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0.02) with the Padova (jGirardi et al. 



derived for three metallicities (Z=0.004, 0.008, 



(120 Id ), blue circles) and with the GALEV 
^Anders fc Fritze-v. Alvensleben ( 2003 ). red triangles) models. The solid curve represents a 
cluster of magnitude rriF555W — 21.5 at different ages at the distance of NGC 4449. For each 
metallicity, the age and mass distributions are also shown (blue solid line for Padova, red 
dotted line for GALEV). 
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fl20r)8h through comparison of the ACS CMD 
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Fig. 19. — Evolution of log(M max .) with equal size (1 dex) log(age) bins for the young (age<l 
Gyr) clusters in NGC 4449, assuming d ifferent metallicities and adopt ing both the GALEV 



( lAnders fc Fritze-v. Alvenslebenl (120031 ) , open circles) and the Padova (IGirardi et al.l (J2010|) , 



full circles) models. Fits to the log(M max ) vs. log(age) for GALEV and Padova are shown 
as dotted and solid lines. Top panels: no internal extinction has been assumed. Central 
panels: we have assumed an internal extinction of E(B — V) = 0.2. Bottom panels: the 
normalized \og(M max /CsFR), where Csfr is proportional to the SFR in the different age 
bins, is plotted. No internal extinction is assumed. See Section [5] for details. 



